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LONG-TERM  GOALS 

Long  term  goals  are  to  obtain  a  thorough  understanding  of  the  behavior  of  the  complete  Stokes  vector 
both  within  the  ocean  and  in  the  atmosphere  as  well  for  both  elastic  and  inelastic  scattering  for  both 
active  and  passive  sources.  Specifically,  we  want  to  see  how  one  can  use  polarization  information  to 
obtain  more  information  about  the  inherent  optical  properties  (IOP’s)  of  both  oceanic  and  littoral  zone 
constituents.  We  also  want  to  explore  the  efficacy  of  polarimetric  time-of-flight  lidar  techniques  in 
determining  salinity  and  speed  of  sound  in  the  ocean  as  well  as  submersible  object  detection.  We  are 
also  working  on  a  new  instrument  design  to  measure  simultaneously  both  the  scattering  and  absorption 
coefficients  of  ocean  water.  We  want  to  continue  our  collaborative  program  in  ocean  and  atmospheric 
polarimetry  with  the  group  in  Minsk  headed  by  Dr.  Eleonora  Zege. 

OBJECTIVES 

It  is  our  major  objective  to  find  new  and  innovative  ways  in  which  polarimetry  can  be  used  to 
determine  not  only  IOP's  but  also  certain  physical  properties  such  as  temperature  and  salinity  of  ocean 
water.  We  also  want  to  further  develop  a  new  method  of  imaging  which  we  would  like  to  label 
"Mueller  matrix  imaging"  (MMI).  With  this  method  we  want  to  determine  how  much  more 
information  we  can  obtain  about  target  features  by  using  elements  other  than  S 1 1  (ordinary  radiance  ) 
of  the  Mueller  matrix  We  will  continue  feasibility  studies  for  a  new  scattering  and  absorption  meter  to 
simultaneously  measure  both  the  scattering  and  absorption  coefficient  for  ocean  water  samples  . 

We  wish  to  continue  our  collaborative  program  with  the  Radiative  transfer  group  in  Minsk,  Belarus 
headed  by  Dr.  Eleonora  Zege.  Our  new  objectives  for  this  program  are  as  follows: 

a)  Rigorously  compare  the  addition  of  a  Lambert  bottom  surface  as  well  as  a  stochastic  interface 
between  the  Multi-component  Approach  (MCA)  and  the  Monte  Carlo  approach  with  full  Mueller 
matrix  formulation. 

b) .  Development  of  the  simplified  equations  for  propagation  of  a  linearly  polarized  beam  in  order  to 
study  the  Stokes  vector  in  the  on-axis  region  for  both  backward  and  forward  directions. 

c)  Development  of  the  code  to  compute  the  Stokes  vector  in  the  on-axis  region  for  a  polarized  narrow 
beam  and  compare  with  experimental  data. 

d)  Extending  the  techniques  we  have  developed  for  detection  and  characterization  of  both  organic  and 
inorganic  atmospheric  aerosols 
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APPROACH 


We  have  developed  some  very  sophisticated  Monte  Carlo  codes  to  handle  the  following  types  of 
problems: 

a)  Complete  Mueller  matrix  calculations  for  a  passive  source  (the  sun)  illuminating  a  plane  parallel 
inhomogeneous  atmosphere-ocean  system  with  a  wind  ruffled  dielectric  interface.  This  program  uses 
a  new  statistical  estimation  technique  that  allows  for  the  multiple  scattering  Mueller  matrix  to  be 
computed  at  precise  angles  without  any  solid  angle  averaging.  We  have  also  developed  a  new 
technique  to  handle  a  Lambertian  bottom  so  that  all  bottom  albedos  can  be  computed  in  a  single  run. 
This  work  is  being  done  by  Hatcher  Tynes,  who  has  recently  completed  his  Ph.D.on  this  project. 

b)  We  have  also  developed  a  new  Monte  Carlo  program  to  calculate  the  multiple  scattered  Mueller 
matrix  from  a  target  embedded  in  a  turbid  medium.  This  program  uses  state-of-the-art  estimation 
techniques  to  improve  the  statistics  over  conventional  "brute  force"  methods.  It  also  utilizes  any 
symmetry  properties  of  the  target.  The  single  scattered  Mueller  matrix  for  the  target  is  a  sine  qua  non 
for  these  calculations.  This  work  is  being  carried  out  by  Deric  Gray  who  is  a  doctoral  student 
working  on  this  project. 

WORK  COMPLETED 

a)  One  phase  of  our  collaborative  work  with  Dr.  Zege’s  group  in  Minsk  has  been  completed  which 
involved  a  comparison  of  her  multi-component  approach  method  with  our  Monte  Carlo  models  for  a 
coupled  atmosphere-ocean  system  with  an  interface  with  no  wind  ruffling  using  a  Henyey  Greenstein 
volume  scattering  function.  We  have  made  the  first  comparisons  of  the  strong  ellipticity  encountered 
underneath  the  interface  of  the  ocean  surface.  The  results  of  this  study  have  recently  appeared  in 
Applied  Optics.1  . 

b)  We  have  completed  a  study  of  the  efficacy  of  using  certain  polarimetric  signatures  to  detect  targets 
at  optical  depths  where  ordinary  radiance  imaging  fails  to  resolve  them.  The  manuscript  for  this 
interesting  work  is  in  preparation  and  will  be  submitted  to  Applied  Optics. 

c)  We  have  also  completed  the  addition  of  the  stochastic  interface  and  Lambertian  bottom  to  our 
effective  Mueller  matrix  Monte  Carlo  Code. 

RESULTS 

We  were  able  to  show  (see  Fig.  1)  that  our  new  statistical  estimation  scheme  could  easily 
accommodate  a  wind  ruffled  interface.  To  check  our  results  we  did  a  comparison  with  the 
HYDROLIGHT  code  of  Curtis  Mobley.  In  Fig.  1  we  show  a  comparison  between  the  two  where 
scattering  in  both  layers  is  totally  conservative  and  is  described  by  the  Rayleigh  scattering  law.  The 
wind  speed  is  5  meters  per  second.  The  atmosphere  has  an  optical  thickness  of  0.01  and  the  ocean  has 
an  optical  thickness  of  1.0.  Sunlight  is  incident  at  the  top  of  the  system,  in  the  direction  of  the  +z  axis, 
which  is  taken  to  be  vertically  down.  Wind  ruffling  is  accounted  for  using  the  statistical  model 
developed  by  Cox  and  Munk.  Below  the  interface,  the  most  noticeable  effect  of  the  wind  ruffling  is  to 
blur  the  boundaries  of  the  so-called  light  cone  (the  cone  that  is  defined  by  the  critical  angle).  As  can  be 
seen  from  the  figure  the  agreement  is  quite  good. 


/  2  2  2 

In  Fig.  2.  We  show  a  plot  of  the  radiance  and  degree  of  polarization  (DOP)=-yQ  +  U  +  V  /I  for  a 
circular  disk  embedded  in  a  scattering  medium.  The  disk  is  at  a  depth  of  4  mean  free  paths  (m.f.p.) 
away  from  a  collimated  source  and  is  a  Lambertian  reflector  with  an  albedo  of  0.025,  but  with  annular 


2 


regions  of  differing  Mueller  matrix  elements  affecting  the  polarization  properties.  1.  Polarization 
preserving  M  =  diag(l, 1,1,1)  2.  Depolarizing  M  =  diag(l, 0,0,0)  3.  Painted  surface  M  = 
diag(l, 0.61, -0.58, -0.51)  The  plot  shows  the  radiance  and  degree  of  polarization  at  a  detector,  located 
at  the  same  position  as  the  source,  for  three  different  input  states  of  polarization  of  the  incident  beam 
relative  to  a  fixed  set  of  axes.  What  is  extremely  interesting  is  that  normal  radiance  detection  would 
not  find  the  target;  however,  by  using  a  single  quantity  such  as  the  degree  of  polarization  (DOP)  we 
can  not  only  "see"  a  hidden  target  but  also  make  out  its  surface  features  when  using  this  quantity  as  an 
imaging  tool  and  more  especially  when  the  incident  beam  is  right  circularly  polarized. 

IMPACT/APPLICATION 

We  are  now  quite  certain  that  Mueller  matrix  imaging,  will  become  a  standard  tool  for  remotely 
sensing  the  atmosphere-ocean  system  as  well  as  an  aid  in  the  detection  of  certain  types  of  cancers. 

We  also  believe  that  it  will  lead  to  better  and  more  comprehensive  target  detection  in  turbid  media. 
What  is  desperately  needed  in  all  MMI  work  is  the  single  scattered  Mueller  matrix  for  the  types  of 
targets  one  is  trying  to  detect. 

Our  Monte  Carlo  programs  for  laser  scattering  from  turbid  media  are  presently  being  used  to  study 
light  scattering  from  living  tissue.  We  are  now  working  in  collaboration  with  Dr.  Ping  Yang  to  expand 
his  FDTD  code  to  handle  biolological  samples  so  we  can  compute  from  first  principles  the  single 
scattered  Mueller  matrix  of  anthrax  spores. 

Another  very  significant  application  of  MMI  is  the  detection  of  biological  aerosols  dispersed  in  the 
atmosphere. 

TRANSITIONS 

Our  Monte  Carlo  passive  source  Stokes  vector  programs  are  being  used  to  check  a  different  method 
for  solving  the  equations  of  transfer  with  full  Stokes  vector  treatment  developed  by  Dr.  Eleonora  Zege 
and  her  group  in  Minsk,  Belarus.  We  are  presently  checking  the  addition  of  a  Lambertian  bottom  as 
well  as  a  stochastic  interface  to  both  our  codes.  Some  of  this  work  was  supported  by  a  CRDF  grant. 
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Fig.  1  The  scalar  radiance  (radiance  calculated  without  including  polarization  effects)  for  an 
atmosphere-ocean  system  with  wind  ruffling  effects  included.  Scattering  in  both  layers  is  totally 
conservative  and  is  described  by  the  Rayleigh  scattering  law.  The  wind  speed  is  5  meters  per  second. 
The  atmosphere  has  an  optical  thickness  of  0.01  and  the  ocean  has  an  optical  thickness  of  1.0. 
Sunlight  is  incident  at  the  top  of  the  system,  in  the  direction  of  the  +z  axis,  which  is  taken  to  be 
vertically  down.  Data  are  compared  between  our  estimation  Monte  Carlo  technique  and  the 
HYDROLIGHT  code  developed  by  Curtis  D.  Mobley.  Wind  ruffling  is  accounted  for  using  the 
statistical  model  developed  by  Cox  and  Munk.  Below  the  interface,  the  most  noticeable  effect  of  the 
wind  ruffling  is  to  blur  the  boundaries  of  the  so-called  light  cone  (the  cone  that  is  defined  by  the 

critical  angle). 
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Fig.2  Plot  of  the  radiance  and  degree  of  polarization  (DOP)for  a  circular  disk  embedded  in  a 
scattering  medium.  The  disk  is  at  a  depth  of  4  mean  free  paths  (m.f.p.)  away  from  a  collimated 
source  and  is  a  Lambertian  reflector  with  an  albedo  of  0.025,  but  with  annular  regions  of  differing 
Mueller  matrix  elements  affecting  the  polarization  properties.  1.  Polarization  preserving  M  = 
diag(l,  1,1,1)  2.  Depolarizing  M  -  diag(l, 0,0,0)  3.  Painted  surface  M  =  diag(l, 0.61, -0.58, -0.51) 
The  plot  shows  the  radiance  and  degree  of  polarization  at  a  detector,  located  at  the  same  position  as 
the  source,  for  three  different  input  states  of  polarization  of  the  incident  beam  relative  to  a  fixed  set 

of  axes. 


RELATED  PROJECTS 

Our  inelastic  scattering  code  is  being  used  in  a  joint  project  with  Dr.  E.  S.  Fry  to  study  the  feasibility 
of  using  Brillouin  scattering  to  measure  the  speed  of  sound  as  a  function  of  depth  in  the  ocean.  This 
project  is  funded  by  the  Texas  Advanced  Technology  Program. 

We  have  also  used  our  time-of-flight  Monte  Carlo  code  to  check  the  validity  of  the  small  angle 
approximation  in  the  solution  of  the  radiative  transfer  equations.  We  have  collaborated  on  some  of 
these  issues  with  Dr.  Walt  McBride  of  NRL  Stennis. 

We  have  also  modified  our  Monte  Carlo  codes  to  emulate  a  new  design  for  a  device  to  simultaneously 
measure  both  the  scattering  and  absorption  coefficients  for  ocean  water  samples. 
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